Imaging of microwave fields using ultracold atoms 
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We report a technique that uses clouds of ultracold atoms as sensitive, tunable, and non-invasive probes for 
microwave field imaging with micrometer spatial resolution. The microwave magnetic field components drive 
Rabi oscillations on atomic hyperfine transitions whose frequency can be tuned with a static magnetic field. 
Readout is accomplished using state-selective absorption imaging. Quantitative data extraction is simple and 
it is possible to reconstruct the distribution of microwave magnetic field amplitudes and phases. While we 
demonstrate 2d imaging, an extension to 3d imaging is straightforward. We use the method to determine the 
microwave near-field distribution around a coplanar waveguide integrated on an atom chip. 



Today, Monolithic Microwave Integrated Circuits 
(MMICs) are of great importance in science and tech- 
nology. In particular, they constitute key building 
blocks of today's communication technology. 1 MMICs 
also serve as main components of superconducting quan- 
tum processors.^ In our group, a simple MMIC structure 
has recently been used as a tool for quantum coherent 
manipulation of ultracold atoms on an atom chipP 

Function and failure analysis is of crucial importance 
for the design of MMICs as well as for simulation 
verification™ External port measurements (e.g. using a 
network analyzer) offer only limited insight. The mi- 
crowave (mw) near-field distribution on the device gives 
much more information, enabling specific improvement. 
Therefore, different methods have been developed to 
measure the spatial distribution of mw near- fieldsPThese 
methods use diverse physical effects to measure the mw 
electric or magnetic field. They have in common that 
they scan the field distribution point- by-point. 

Here we propose and experimentally demonstrate a 
highly parallel method that allows for non-invasive and 
complete (amplitudes and phases) imaging of the mw 
magnetic field distribution using clouds of ultracold 
atomsP In this method, the mw magnetic field drives 
resonant Rabi oscillations^ between two atomic hyperfine 
levels that can be detected using state-selective absorp- 
tion imagingP The method offers fim spatial resolution 
and a mw magnetic field sensitivity in the 10 -8 T range 
at frequencies of a few GHz. It is a frequency-domain, 
single-shot technique to measure a 2d field distribution. 
The method can be extended to measure 3d distributions 
slice by slice. Data extraction is simple, it offers a high 
dynamic range, and it is intrinsically calibrated since only 
well-known atomic properties enter in the analysis. 

For the proof-of-principle experimentpresented here, 
we use our atom chip setuppl see Fig. [TJ The chip has 
integrated mw coplanar waveguide structures that were 
designed for quantum manipulation of ultracold atoms P 
To demonstrate our method, we analyze the mw mag- 
netic field near this structure. In general, the device to 
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be tested does not have to be integrated on an atom chip. 
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FIG. 1. Schematic of the experiment, (a) Ground state hy- 
perfine levels of 87 Rb atoms in a static magnetic field. Ini- 
tially, the atoms are trapped in state \F, tuf) = |1, —1). The 
three relevant transitions |1, — 1) — > \2,m2), (m2 = —2,-1,0) 
are indicated. The corresponding transition frequencies w 7 , 
(7 = — ,7r, +) are split by ui-l due to the Zeeman effect. The 
Rabi frequencies fi 7 are indicated, (b) The atom chip. Inset: 
Atom cloud near the coplanar waveguide (CPW) structure 
whose mw magnetic field is examined. The three inner wires 
constitute the CPW. (c) Experimental sequence. Left: The 
trap is switched off and the atom cloud expands. Right: A mw 
pulse is applied to the CPW, resonant with one of the transi- 
tions uj~f. Its magnetic field of amplitude B(r) drives Rabi 
oscillations with position-dependent fi 7 (r) between |1, — 1) 
(red) and the corresponding state 2,7712) (blue). The re- 
sulting atomic density distribution ni(r) (nz(r)) in F = 1 
(F = 2) is detected. From n;(r), we reconstruct f2 7 (r). Sev- 
eral such measurements on the three transitions are combined 
to reconstruct B(r). 

Our method works as follows. On the room- 
temperature atom chip, we prepare clouds of Rb atoms 
at a temperature of T ~ 5 /zK using laser and evapo- 
rative cooling techniques!^ The magnetically trapped 
atoms are initially in the ground state hyperfine sub- 
level \F, rap) — 1 1,-1), see Fig. [iJl The trap is moved 
close to the mw structure to be characterized, where it is 
switched off and the atoms are released to free fall. Dur- 
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ing a hold-off time eft ho, the cloud drops due to gravity 
and expands due to its thermal velocity spread, filling 
the region to be imaged (Fig. lb+c). We maintain a 
homogeneous static magnetic field Bo of order 10~ 4 T. 
It provides the quantization axis and splits the frequen- 
cies u>y, (7 = — ,7r, +) of the three hyperfine transitions 
|1,-1) -> |2,m 2 ), (m 2 = -2,-1,0) by cj l = fi B B /2H 
due to the Zeeman effect (cf. Fig. la). A mw signal on 
the MMIC is subsequently switched on for a duration 
di mw (typically some tens of fis). We select one of the 
transitions by setting the mw frequency u> = cj~. The 
mw magnetic field S§{r,t) = \ [B(r)e" lLJt + B*(r)e iwt ] 
couples to the atomic magnetic moment and drives Rabi 
oscillations of frequency ^i 7 (r) on the resonant transition. 
For the three transitions of interest, 
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Here, -B,r and are the real-valued amplitude and 
phase of the component of B parallel to Bo, and B + , (j> + 
(-B_,0_) are the corresponding quantities for the right 
(left) handed circular polarization component in the 
plane perpendicular to Bo, see Supplementary Informa- 
tion. After the mw pulse, a spatial pattern of atomic 
populations in F = 1 and F = 2 results, see Fig. [I]:. The 
probability to detect an atom in F = 2 is 
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Here, ni(r) (n 2 (r)) is the density of atoms in F = 1 
(F = 2), which can be measured using state-selective ab- 
sorption imaging.® From p 2 (r) we can reconstruct |f2 7 (r)| 
and thus the spatial distribution of the resonant mw po- 
larization component £? 7 (r), as shown below. 

For a reconstruction of B(r), we measure B ff (r), 
B + (r), and B_(r) with B pointing along the x, y, and 
z axis. In each of the nine measurements (cf. Fig. [2]), 
we tune ui into resonance with the desired transition w 7 . 
Alternatively, to test a device at a given frequency u, 
one can achieve resonance by Zeeman-tuning of cj 7 via 
Bq . The measurements allow a reconstruction of the am- 
plitudes of all three cartesian components of B as well 
as their relative phases. The measurements of B„ for the 
three orientations of B directly yield the amplitudes B x , 
By, and B z . The relative phases can be reconstructed 
from the other measurements, see Supplementary Infor- 
mation. It is also possible to measure the spatial depen- 
dence of the global phase of B using an interferomctric 
method, see Supplementary Information. 

Absorption imaging integrates over the direction of 
propagation of the imaging laser beam. Structures to 
be characterized should therefore have a characteristic 
length scale along the beam larger than the size of the 
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P m „=120 mW, dt m „=20 ms P m „=120 rnW, dt m „=20 jjs P m „=120 mW, dt m „=40 [js 

FIG. 2. Imaging of mw magnetic field components near the 
CPW of Fig. [l] The images show the measured probability 
Pa(r) to find an atom in F — 2 after applying the mw pulse. 
Columns correspond to measurements on the three different 
transitions o; 7 , rows to three different orientations of Bo. The 
imaging beam is reflected from the chip surface at an angle of 
2°. As a result, on each picture, a direct image and its reflec- 
tion on the surface are visible. The dashed line separates the 
two. No atoms are visible in the very center of the image be- 
cause the CPW structures distort the imaging beam. The mw 
power launched into the CPW, P mw , and the mw pulse dura- 
tion dimw are indicated, dtho varies between f - 2 ms. The 
noise in the image periphery corresponds to regions without 
atoms. Images are averaged over several experimental runs 
(15 to 130). 



atom cloud. The cloud size can be adjusted through the 
magnetic trap frequencies, T, and e?iho- 

We experimentally demonstrate our method by mea- 
suring the mw magnetic field distribution near the copla- 
nar waveguide (CPW) on our atom chip. The measure- 
ment is performed at a position where the CPW is trans- 
lationally invariant along the imaging beam. An overview 
of the data is shown in Fig. [2j What appears in the im- 
ages are essentially isopotential lines of the mw magnetic 
field components. Qualitative conclusions can be drawn 
directly by looking at the images - e.g. the left /right 
asymmetry visible in Fig. [2] reveals that there is an asym- 
metry in the mw currents on the CPW wires. This con- 
firms independent findings in our previous paper P 

It is possible to automatically extract B 1 (r) from a 
single image. From p 2 (r) we can calculate |fi 7 (r)| using 
Eq. Q up to an offset of n7r/o?t mw , where n is an in- 
teger. The offset for each point can be calculated by a 
ray-tracing method, where rays are sent from the image 
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FIG. 3. Measured mw magnetic field component B x (r) near 
the CPW and comparison with a simulation. The measure- 
ment is performed with ui = cj^ and Bo along x, and the 
movie method (see text) is used to extract B n (r) = B x (r). 
(a) Rabi oscillations at two exemplary pixels of the image, 
recorded by varying P mw at fixed dt mv = 20 /j,s. The sinu- 
soidal fits used to determine |fi,r| and thus = B x as a 
function of \Afmw are shown. The observed decay of the os- 
cillations is due to mw field gradients across the pixel, (b) 
Image of B x (r) at P mw = 120 mW as obtained from the data, 
(c) Corresponding quasi-static simulation of B x (r)W We find 
best agreement with (b) if we allow for a 10% asymmetry be- 
tween the currents on the two CPW ground wires and assume 
induced currents in the two wires next to the CPW grounds 
with an amplitude of 2% of the signal conductor current. 



periphery (where \tt 7 \dt mv/ -C 1) through the desired 
point to the center of the mw structure (where |fi 7 | is 
maximal), n is given by the sum of the number of min- 
ima and maxima of P2 (r) encountered on the ray. 

Alternatively, we take series of k image frames, scan- 
ning either the mw power P mw or dt mw (the movie 
method), k depends on the desired dynamic range, 
but can be as low as 10. The time to record one 
frame is 14 s, but could be reduced to < 3 sP^ For 
each image pixel, we thus obtain a sequence of k dat- 
apoints showing Rabi oscillations, see Fig. [3^,. We fit 

) | dtmw] to the data, where 
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and a is the fit parameter. From 




the fit, we determine |f2 7 | and thus, via Eqs. 
B 1 at this pixel for a given P mw . As an example, Fig.|3|: 
shows an image of the cartesian mw field component 
B x (r) near our CPW reconstructed in this way. 
shows a corresponding simulation. Comparing data and 
simulation, we obtain information about the current dis- 
tribution on the CPW. 

The single-shot mw field sensitivity of our method is 
mainly determined by the interaction time dt mw . There is 
a trade-off between sensitivity and effective spatial res- 
olution s c g. Longer dt mw yields higher sensitivity, but 
at the same time the image blurs due to the movement 
of the atoms. For our parameters (optical resolution 
s opt = 4 /im, dt mv/ = 80 fis, T — 5 /iK) and dt\ w = 0, we 
obtain s e g = 8.2 /im and a mw magnetic field sensitivity 
of 2 x 10 -8 T (Supplementary Information). 

In a variant of the presented method, trapped atoms 
could be used as a scanning mw field sensor. For trapped 
atoms, dt mw can be much longer without increasing s e g , 
thereby improving the field sensitivity. To record an im- 



age, the trap position has to be scanned from shot to 
shot of the experiment (Supplementary Information). 

Our method is frequency selective. Individual com- 
ponents of a multi-tone signal could be resolved. The 
transition frequencies w 7 can be adjusted via Bq. For 
Bq up to 0.5 T, which is a realistic effort e.g. for proto- 
type testing, transition frequencies of 2.5 - 14 GHz are 
accessible with 87 Rb (Supplementary Information). Note 
that for Bo > 0.1 T we start entering the Paschen-Back 
regime, where the matrix elements of Eqs. ([!]) - ^ change 
and the theory has to be modified. Using other atomic 
species, different frequency ranges are accessible. Alter- 
natively, a two-photon transition could be used for imag- 
ing, where two mw fields or a mw and a radio frequency 
field are applied to the atoms. The first field (frequency 
lui ) is applied externally with known spatial distribution, 
while the other field (frequency w 2 ) is to be imaged. Res- 



onant Rabi oscillations occur for uj\ + uj-2 



It is also 



possible to image off-resonant mw fields and other dif- 
ferential potentials using a Ramsey interferometry tech- 
nique, see 8 ^ and Supplementary Information. 

Our technique can be extended to measure 3d distri- 
butions of B(r) slice by slice, either by using a gradient 
of Bq such that only a slice of atoms is resonant with 
uj, or by using a light sheet detection technique,^ where 
slices perpendicular to the camera line of sight are im- 
aged. Using similar techniques, it is also possible to shape 
the atomic cloud to prepare a thin sheet of atoms. 

The method presented here seems promising for ap- 
plications like prototype characterization. It allows 
for highly sensitive, parallel, high-resolution, and non- 
invasive imaging of the complete mw magnetic field dis- 
tribution around an MMIC. Compact and portable sys- 
tems for the preparation of ultracold atoms have been 
built P and key components of such systems are com- 
mercially available. 
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IMAGING OF MICROWAVE FIELDS USING ULTRACOLD ATOMS: SUPPLEMENTARY INFORMATION 



A. Rabi Frequencies 

We derive the Rabi frequencies for the resonant coupling of ground state hyperfine levels of 87 Rb with a microwave 
(mw) field. In the following, we consider an atom in a weak static magnetic field Bo, so that the Zeeman splitting 
k>£ = is small compared to the zero-field splitting of ujq ~ 2tt x 6.8 GHz between the two hyperfine states 

F = 1 and F = 2 of the 5 2 Si/2 electronic ground state of 87 Rb. The atom is initially prepared in the hyperfine 
sublevel \F,mp) = |1, — 1), and the microwave frequency w is resonant with one of the transitions | X, — 1) O |2, m 2 ), 
(to2 = —2,-1,0) connecting to this level (see Fig. 1 of the main paper). 

The real-valued mw magnetic field at position r = (x, y, z) in the fixed cartesian laboratory coordinate system is 
38(r, t) = \ [B(r)e _lwt + B*(r)e lLUt ] with the complex phasor 

/ B x (r)e-^^ 
B(r) ee S y (r)e-^W 
\ B,(r)e-*W 

Here, we have chosen £?j(r),0j(r) £ K>o, {i — x,y,z). In the following, we consider a fixed position in space and 
suppress the dependence of Bi(r) and 4>i(r) on r to simplify notation. In our experiment, we apply a homogeneous 
static magnetic field Bo along several directions in order to be able to reconstruct all components of the mw magnetic 
field. For a given B , we choose a new cartesian coordinate system (x',y',z') with the z'-axis pointing along B , 
which defines the quantization axis for the atomic states \F, mp). In this new coordinate system, the mw magnetic 
field phasor is given by 

B = | By.e-^v 
B z ,e-^ 

The mw magnetic field couples to the magnetic moment of the electron spin of the atom. The coupling to the nuclear 
magnetic moment is neglected, because it is three orders of magnitude smaller than the electron magnetic moment. 
The Rabi frequency on the hyperfine transition |l,mi) f-> |2, 771.2) is given by 

^ = 2 ^<2,m 2 |B.J|l,m 1 ), (5) 

with J — (J x >, J y ' , J z ') the electron spin operator. Using J± — J x r ± iJ y i we can write 

B J = B x ,e'^' J x > + B y ,e-^y' J y , + B z ,e~^*' J z , (6) 

= \ (B x ,e~ 1 ^' -iBye'^y') J+ + \ (B x ,e'^' + %B v ,e^y') J- + B^eT** J z ~. (7) 
Evaluating the matrix elements for the three transitions connecting to 1 1, — l)p we obtain the Rabi frequencies: 

q_ = n a,-a = 2^b (2j _ 2| 1 ( Bx , e -i^, + iByie -i^) _i) = _ e -#- . ^3 . ^ B _^ (8) 

iK^nlzl = 2 ^(2,-i|^ e - i ^'j^|i,-i) = -e-^- ■]jl- > Y B «, (9) 

n+^tf^ = ^(2,0|i (B^e-** -iB y ,e-^') J+|l,-l) = e"^ ■ ^B + , (10) 
where we have used the definitions 

B_e-**- := § (B^e-'^' + iB^e-***') , (11) 

B T e-*+* ^ Bje-***' , (12) 

B+e-^+ := § (B^e-^' - iBye"^') , (13) 

with B 1 ,4> 1 £ K>oj (7 = ~ 7 71 ", +)• We note that f2„. is proportional to the projection of B onto B , while f2+(_) is 
proportional to the right (left) handed circular polarization component in the plane perpendicular to Bo- 

In the experiment, we choose a sufficiently strong static field so that lul 3> Q 7 . Furthermore, we choose the 
microwave frequency resonant with one of the transitions, u> — cj 7 . In this way, Rabi oscillations are induced only on 
the resonant transition in a given run of the experiment, which allows us to selectively image the individual microwave 
magnetic field components _B 7 . 
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B. Reconstruction of the microwave magnetic field 



The amplitudes B x , B y , and B z of the cartesian components of B in laboratory coordinates can be easily determined 
by measuring |f2„.| with the quantization axis Bo/Bo pointing along x, y, and z, respectively. The extraction of the 
field components from absorption images is described in the main text. In the following, the upper index indicates 
the direction of the quantization axis in laboratory coordinates, e.g. Q v _ (B v _) means f2_ (£>-) for Bo pointing along 
the y-axis. 

To reconstruct the relative phases (4> y — <f> x ) and (</> z — <j) x ) between the cartesian components of B, we also measure 
the amplitudes of the circularly polarized components B^, B^_, B+, B v _ 7 B+, and B z _- Having measured these 
components, we can reconstruct the relative phases according to the following recipe. 



1. B along x 

We choose a new coordinate system with z' along x, resulting from the following coordinate transformation: 



x' 


= —z 


y' 


= y 


z' 


— X 



In this rotated coordinate system, the mw magnetic field phasor reads 

-B z e 
Bye'^y 
y B x e- i4 >* 

From this we obtain: 



B 




B x > 


= B z 


4>x> = 4>z + n 


By, 


= By 


4>y> = <t>y 


B z , 


= B x 


4>z> = 4*x 



Using Eqs. (11) and (13), we calculate 

B\ - B 2 _ = —B x ,B y i sin(<^, - 
By insertion of the coordinate transformation and using Eqs. 

h 



^ B B z B y V3 



I SI 



x |2 



n 



x |2 



( 10 ), we obtain 
1 



flByB Z 



((B*_f-(Blf) 



(14) 



(15) 



2. B along y 

A similar calculation as before yields 

sm((j) x - (f) g ) 



^i B B x B z \Z 



in* 



8|fi 



y i2 



((Blf - (Blf) 



(16) 



3. Bo along z 

In this case, we obtain 



sin ( 



4^ B B x B y V3 



10? 



10? 



(17) 



All quantities on the right hand sides of Eqs. (15) - (17) can be measured. From Eqs. (15) - (17), the relative 
phases {<f> y — 4> x ) and (<fi z — 4> x ) can be determined. The solution is unique except for the very degenerate case where 
sin (4> x — <j) z ) = sin [<f> y — <f> x ) — sin (4> z — 4> y ) — 0. In this case, there are 4 solutions which cannot be distinguished. 
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C. Reconstruction of the absolute microwave phase 

While Sec. [B] describes the reconstruction of the relative phases (4> y — 4> x ) and (4> z — 4> x ) between the cartesian 
components of B, it is also possible to reconstruct the spatial dependence of the global phase of B. The procedure 
uses two mw pulses. During the whole sequence, B stays the same. In the following, we assume B is pointing 
along the x-axis, so that <fi x (r) is measured. The other two phases, 4> y (r) and <^ z (r), can then be determined from the 
already known relative phases. 

After releasing the atoms from the trap, they are prepared in an equal superposition of states |1,— 1) and |2, — 1) 
by application of a ^-pulse at frequency u) = This mw pulse is applied from a well-characterized source, so that 
it has negligible (or at least known) intensity gradients and negligible (or known) phase gradients across the atomic 
cloud. This can be achieved by using an external mw hornP The duration of the pulse is di mw i = girf i| > wnere 

f^i = e 1 ^' 1 is the Rabi frequency of the pulse. The state after this preparation pulse is (in the rotating wave 

approximation^ 

IV>i) = ^(|l,-l)+*e-^|2,-l}). (18) 

Immediately after the end of this preparation pulse, the mw in the MMIC to be characterized is pulsed on at frequeny 
<^mw,2 = <^tt for a duration dt mw ,2- The Rabi frequency and phase of this second mw pulse are denoted by f2 ff ,2(r) 
and <^7r,2(r), respectively. After the second pulse, the state of an atom at position r is 



l^(r)) = 4 



2 ) \ 2 

ie #.,i cos f l^2(r)|^ mw , 2 \ +ie -i^,= si _ f K,2(r)|dt ro w,2 



1-1) 

2,-1). 



s/2 

+ 73 

The probability P2 (r) of finding an atom at position r in state F — 2 is given by 

P2(r) = - + -sin(|fi 7r ^ 2 (i")|^mw,2) • COS (0^,1 - <pw,2( T )) ( 19 ) 

To calculate Wj 2( r ), the quantities |f2 w 2( r )| an d 4>n,i have to be known. |f2„. 2 (r)l can be measured as described in 
the main text. If Bo is pointing along the x-axis, then <j) x (r) = <f> Wt 2(r)- 

The calculation above assumes that there is zero delay between the end of the first preparation pulse and the second 
pulse in the MMIC. A similar calculation is also possible for non-zero delay between the two pulses. 



D. Sensitivity & spatial resolution 

In this section we estimate the maximum sensitivity of our technique for our set of experimental parameters. As 
described in the main text, the mw magnetic field sensitivity is mainly determined by the interaction time dt mw of 
the atoms with the mw pulse on the MMIC. A longer interaction time dt mw leads to a higher mw magnetic field 
sensitivity, because then a weaker mw field can already drive a substantial fraction of a Rabi cycle. However, at the 
same time the effective spatial resolution s c g = 2a Q g decreases as the image blurs due to the movement of the atoms 
during dt mw . In the following, the symbol a always refers to an r.m.s. width. a c g is determined by the average moving 
distance er mw of the atoms during the mw pulse, by the optical resolution of the imaging system s opt = 2a op t , and by 
the movement of the atoms during the imaging laser pulse, where thermal motion (<r tm ) and diffusive motion due to 
photon scattering (<7 ps ) contribute. In the following, we will calculate tr mw , <r tm , a ps , and <7 e ff. 



1. Movement of atoms during the microwave pulse - <r mw 

A free-falling atom in a cloud at temperature T has a mean thermal velocity perpendicular to the line of sight of 
v t h = y 2 *^ r ^ After releasing the atom from the trap and waiting for dt\ w , the atom has furthermore acquired a 
velocity v g — g ■ dt\ lo along the direction of gravity. During the interaction with the mw for a time dt mw , the atom 
moves ballistically by an average distance 



1 / 2k T 

fmw = 9 ■ dt ho ■ dt mw + -g ■ dt w + \ dt mw . (20) 

z \ m 



For T = 5 /iK, dt\ la = 0, and dt mw = 80 /is we obtain er mw = 2.5 /im. For short times, the last term in the equation 
above dominates (as it is the case for our parameters). The displacement of the atoms during the mw pulse can then 
be approximated by a Gaussian function / mw of r.m.s. width <r mw . 



2. Movement of atoms during the imaging pulse - a t: 



During the imaging pulse of duration dt- lla , the atoms move ballistically by an average distance cr tm given by 



1. _u2 , WET 



0tm = g(dt ho + dt mw )dt im + -g ■ dt im + \ • dt im (21) 

2 \ m 

due to gravity and thermal motion. For dt[ m = 40 [is we get Otm = 1-3 /im. Again, for short times the atomic density 
distribution after the imaging pulse can be approximated by a Gaussian / tm with r.m.s. width <7t m . 



3. Diffusive movement of atoms due to photon scattering - a ps 



During the imaging laser pulse of duration dt lra , the atoms randomly scatter photons. The associated momentum 
recoils lead to a diffusive motion of the atoms, which leads at the end of the pulse to an average displacement 
perpendicular to the line of sight oP 




where v ICC = hk/m = 5.9mm/s is the atomic recoil velocity for 87 Rb, and N p the number of scattered photons, with 
Np = (r/2) dti m s/(l + s) and the natural linewidth r = 2tt x 6.1 MHz. For our experimental parameters (saturation 
parameter s = I/I sa ,t = 1, dt lm = 40 fis) we get <7 ps = 2.2 /im. 



4. Effective spatial resolution s e g = 2a c g 



The effective resolution can approximately be calculated by the convolution f c g = / opt * [/ ps * (/tm * /mw)]> where 
/opt approximates the point spread function of the imaging system by a Gaussian / opt with a op t = 2 /im. As a result, 
we get a e g = 4.1 /im. We take as the effective resolution for our parameters s c g — 2a e g = 8.2 /im. 



5. Imaging noise 



The noise on the absorption images is important for the sensitivity of our technique. It determines the minimum 
number of atoms A^min that has to be transferred by the mw into the F = 2 manifold during dt mw in order for the mw 
field to be detectable. We currently use an Andor iKon-M camera with a quantum efficiency of 90% for absorption 
imaging^. The optical resolution of our imaging system is s pt = 4 ^m, the imaging pulse duration dt- lm = 40 [is, and 
we are imaging at saturation intensity (s = I/I aa ,t — 1)- With these parameters, we calculate an uncertainty in the 
number of atoms detected in an area A e g = 7rcr^ ff of 1.4 atoms r.m.s. We measure a value of <Tjv, pS n = 2.0 atoms. The 
difference can be explained by interference fringes on the image. This additional noise could certainly be decreased 
further. 

Quantum projection noise due to the probabilistic nature of the measurement process is an additional contribution 
of noise on the images. The measurement process projects the atomic superposition state onto the F = 1 and F = 2 
states, resulting in a number of atoms of N\ and N% in the two states, respectively. Even if the total number of atoms 
N = Ni + N2 is the same in each shot, N\ and N2 will show (anticorrelated) fluctuations. This projection noise has an 
r.m.s. amplitude of cr/vi = <7/v2 = y^V • P2 ■ (1 — P2), where cr^i denotes noise in Ni and P2 = N2/N. The total noise 

on N2 is thus u jv,tot — psn + a N2 ■ We find that in order to obtain a signal-to- noise-ratio SNR = A^/crjv.tot > 1, 
we have to have N2 > i^min = 3. 
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6. Microwave field sensitivity 

In the experiment described in the main paper, we trap about N = 9 x 10 3 atoms in the magnetic trap. The 
trapping frequencies are uj x = 2n x 27 Hz and uj y ~ lj z — 2ix x 680 Hz. We calculate the average atomic density in the 

trap to n = 2.2 x 10 11 /cm 3 PThe trapped cloud has a 1/e-radius of p = \J 2k ^ F ^- = 7.2 /im along the y-axis, which 
is the direction of the imaging beam. If we image the atoms with dtho = 0, we have about N — 170 atoms inside a 
cylinder of radius a c g and height 2p. The mw magnetic field which transfers on average A^min = 3 atoms to F = 2 
is obtained by requiring that 

= N 2imin . (22) 
For the three transitions uj = w_, w = w„ and U) = w+, this is equivalent to 

2 

= 3, (23) 



N 2 = N sin 2 



|f2 7 |dt r 



N F=2 = N sin 2 



N 



N, 



F=2 



TV sin 2 



Ms 

h 



-B« dt 
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N 



fJ-B 

h 



-Bjr dt 

n 



3. 



(24) 



N 



F=2 



TV sin 2 




N 




(25) 



= 80 fxs, we get B_ = 2.2 • 10~ 8 T, B^ = 4.4 • 10~ 8 T, and B + = 5.4 • 10~ 8 T. 

Therefore, the absolute mw field resolution 



Solving Eq. (22) 



Solving the above equations with dt 

Note that the projection noise ctjv2 slowly increases with increasing N2 
of our method decreases with increasing values of B. 

The consideration above relies on the assumption that we can achieve perfect resonance uj = oj 7 
for A^min, we obtain |J7 7 |/27r = 0.53 kHz for 7V2 = 3. A change in w 7 /27r of 0.53 kHz corresponds to a magnetic 
field instability of 2.5 x 10~ 8 T for 0_, 3.8 x 10~ 8 T for ^ and 7.6 x 10~ 8 T for fl + . Inside the magnetic shielding 
surrounding our experiment, we achieve a stability of Bq of 2 x 10~ 8 T r.m.s., which could certainly be improved such 
that the effect can be neglected. 

The sensitivity of our field imaging technique can be increased by using colder or denser clouds. Suitable techniques 
to reduce the temperature further are adiabatic relaxation of the trap or further forced evaporative cooling. 



E. Measurement of the microwave field with trapped atoms 



There are further operation modes of our field imaging method. Instead of releasing the atoms from the trap, it 
is also possible to prepare a Bose-Einstein condensate (BECjf^ in a trap (either magnetic or optical) or in an array 
of traps (e.g. an optical lattice) and scan its position spatially from shot to shot. The position can be scanned on 
a /im length scale in all three dimensions. A typical stepsize is given by twice the Thomas-Fermi radius ttf of the 
BECPlf the BEC is held in a magnetic trap, transitions between |1, — 1) and |2,m2) are detected via loss of atoms, 
because |2, m 2 ) with ni2 = —2, —1,0 are magnetically non-trappable states. If the BEC(s) are held in an optical trap 
or optical lattice, then transitions to the F = 2 manifold can be detected via state-selective imaging. 



F. Tunability of transition frequencies oj , cj n and lo + 
1. 87 Rubidium atoms 



The transition frequencies w_, lu^, and ui + between the initial state |1, — 1) and the target states |2,— 2), |2, — 1), 
and 1 2 , 0) , respectively, can be adjusted by changing the static magnetic field B^. For small magnetic fields (i.e. 
ujl <^ 2tt x 6.8 GHz) the transition frequencies are approximately given by w_ — ujq — 3w^, uj^ — uj — 2ll>l and 
ui + = ojq — ujl, where ujq — 2tt x 6834.682610 MHz and lul = ^bBq/2H. 
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For larger values of Bo a more accurate calculation is required. The 5 2 Si/2 ground state hyperfine structure of 
87 Rb is described by the HamiltoniarP 



H = A Ms I ■ 3 + LiBB Q (gjJ z + gil z ). 



(26) 



with I the nuclear spin operator and J the operator for the total electronic spin, Am s = Hujq/2, gj = —0.000995141, 
and gj = 2.002331.-J Even though |y « 5 • 10~ 4 , the coupling of I to the magnetic field becomes important for large 
values of Bq. An analytical formula exists for the ground state manifold of a D transition (as it is the case for our 
states of interest), the Breit-Rabi formula^: 



E 



A 



hl's 



F,m j 



(2J + 1) 



gi\i B m F B ± A his 1 



4m^ 2 
21 + 1 M 



1/2 



(27) 



where (3 — ty9J 2A^ bB ° 1 m F — mi±mj (the ± sign is the same as in eq. 27). The resulting transition frequencies cj_, 
llv, and u> + as a function of B are illustrated in Fig. |ij For 87 Rb, the transition frequencies w 7 /27r can be tuned over 




2 3 4 

Bo[10" 1 T] 

FIG. 4. Transition frequencies ui n and ui+ for 87 Rb as a function of -Bo- The individual transition frequencies can be tuned 
by up to 10 GHz with magnetic fields of up to 0.5 T. 



a range of more than 10 GHz using technically feasible magnetic fields. Note that for B > 0.1 T we start entering 
the Paschen-Back regime, where the matrix elements of Eqs. (p| - ( 10 1 change and the theory has to be modified. 



2. Other atomic species 



By using atomic species other than 87 Rb, different frequency ranges become accessible, e.g. 9.2 GHz for Cs or 
1.7 GHz for Na. 



G. Ramsey interferometry and off-resonant probing 



Instead of having lo resonant with a hyperfine transition frequency, it is also possible to probe an off-resonant mw 
or light field (or anything else that causes a differential energy shift between the involved hyperfine levels) using a 
scheme based on Ramsey interferometryP In Ramsey interferometry, the interaction with the mw field of duration 
Tj nt is enclosed by two Impulses. The first pulse prepares the atoms in an equal superposition of two hyperfine states 
such as 

hh>(* = 0) = \ (1 1,-1) +ie- i *™|2,m 2 )) . (28) 



11 



During time Tint both states accumulate a differential phase shift Acj) = J Q mt ^jj^dt. E^is is the differential potential 
between the two hyperfine levels involved and is in general state-dependent. The state after this interaction is 

l^i>(* = T int ) = ^ (|1, -1) +ie- iA ^»|2,m 2 )) . (29) 

After applying the second ^-pulse, the state is 

|V>2> = \ [(1 - e-^) |1, -1) + ie-*™ {e- 1 ^ + l) \2, ma)] , (30) 

where mw is the phase of the resonant mw in a frame rotating at the atomic transition frequency. The probabilities 
to detect an atom in state |1, — 1) and |2, ra-i) are then given by 

Pi(r) = i(l-cos(A0)), (31) 
P2(r) = l-pi(r) = i(l+co8(A0)). (32) 

By measuring the relative populations pi(r) and P2(r) after the second ^ pulse, it is possible to determine the value 
of A(f> and thereby -Ediff- Details on the effect of an off-resonant mw field are discussed irP. 
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